An area law for the entropy of low-energy states 
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It is often observed in the ground state of quantum lattice systems with local interactions that the 
entropy of a large region is proportional to its surface area. In some cases, this area law is corrected 
with a logarithmic factor. This contrasts with the fact that in almost all states of the Hilbert space, 
the entropy of a region is proportional to its volume. This paper shows that low-energy states have 
(at most) an area law with the logarithmic correction, provided two conditions hold: (i) the state 
has sufficient decay of correlations, (ii) the number of eigenstates with vanishing energy-density is 
not exponential in the volume. These two conditions are satisfied by many relevant systems. The 
central idea of the argument is that energy fluctuations inside a region can be observed by measuring 
the exterior and a superficial shell of the region. 
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I. INTRODUCTION 

Entropy quantifies the uncertainty about the state of 
a physical system. A bipartite system in a pure state 
has zero entropy, but the reduced state of one subsystem 
may have positive entropy. This is due to quantum cor- 
relations between the two subsystems, the entanglement. 
In fact, this entropy quantifies the entanglement in the 
sense of quantum information theory [l[ . 

In classical physics, the entropy of a region inside a 
spatially-extended system at finite temperature is pro- 
portional to the volume of the region — entropy is an ex- 
tensive quantity. At zero temperature, it is small and 
independent of the region. In quantum physics, at finite 
temperature, the entropy of a region is also proportional 
to the volume. But it has been observed in several mod- 
els that, at zero temperature, the entropy of a region is 
proportional to its surface area (2l4To|. In some mod- 
els of critical free fermions the entropy scales as the area 
times the logarithm of the volume [ill, E3 • This has been 
presented as a violation of the area law, although the di- 
mensionality of the scaling is still that of the area. A 
celebrated proof shows that any one-dimensional system 
with finite-range interactions and an energy gap above 
the ground state obeys a strict area law (l3| . 

The original motivation for this problem is the anal- 
ogy with black-hole physics, where the thermodynamic 
entropy is proportional to the surface area of the event 
horizon @, 0, [l4| ■ The second motivation is to guide the 
development of efficient methods for simulating quantum 
systems with classical computers. The number of param- 
eters needed for specifying an arbitrary pure state of an 
TV-partite system is exponential in N. If the state is not 
entangled, the number of parameters is proportional to 
N. Hence, there seems to be a correspondence between 
entanglement and complexity. In one spatial dimension, 
the relation between entropy and the complexity of sim- 
ulating a system is well understood [B|, [HI, [l6| . The third 
motivation is to understand the kind of states that arise 
in quantum many-body systems with strong interactions. 
Almost all states in the Hilbert space obey a volume law 
for the entropy [TtJ ■ Hence, area laws tell a lot about the 
multipartite entanglement structure. At a finer level, the 



specific form of an area law tells additional information 
about the system: the logarithmic correction is a signa- 
ture of criticality [1, 0, 

H EH, E3; and the appearance of a 
negative constant is a signature of topological order [l8[ . 
For further overview of the topics around area laws see 
the review article cited [l9| . 



II. RESULTS AND SUMMARY 

Consider an arbitrary hamiltonian H with finite-range 
interactions in an s-dimensional lattice. The eigenstates 
have a well-defined global energy, but inside a region X 
of the lattice the energy may fluctuate. (The nomen- 
clature of FIG. 1 is followed.) In Section III it is proven 
that these fluctuations can be observed by measuring the 
exterior of the region and a superficial shell inside the re- 
gion, that is XUS. In Section IV a condition is imposed 
to the ground state: if the operator X has support on the 
region TZ which is separated from the support of the op- 
erator Y by a distance I, then the connected correlation 
function decays at least as 



\(XY}-(X)(Y)\<(l-£ln\K\)- 



(1) 



where £ is a constant. This implies that energy fluctua- 
tions inside the region X cannot be observed in its bulk, 
namely TZ. This provides a characterization for the ap- 
proximate support of the global ground state inside the 
region TZ. In Section V a condition on the density of 
states is assumed: if Hx is the subhamiltonian with all 
terms of H whose support is fully contained in X, then 
the number of eigenvalues lower than e is bounded by 



Ov(e) < ( r \X\)^ e - e ^ +r '\ dx \ 



(2) 



where eo is the lowest eigenvalue and r, 7, 77 some con- 
stants independent of X. This condition is only assumed 
for e ~ \dX\. This implies an upper-bound on the dimen- 
sion of the above-defined support subspace. This is used 
to bound the Von Newmann entropy for the reduction of 
the global ground state in the region 1Z 

S{pn) = tr(-p n hxpn) < const \&R\ ln\K\ . (3) 
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where the index n = 0,1, 
increasing oreder E n < E n . 



. . labels the eigenvalues in 



Note that any hamiltonian with finite-range interac- 
tions in a sufficiently regular lattice can be brought to 
the form of H, by coarse-graining the lattice. Quantum 
field theories with local interactions can also be brought 
to the form of H by lattice regularization. In the case 
of bosons, a truncation in the local degrees of freedom 
is needed. In the case of fermions, a multi-dimensional 
Jordan- Wigner transformation [2(| is needed. 



The Lieb-Robinson Bound 



FIG. 1: (Color online) 1Z is the chosen region where the en- 
tropy is estimated; the sites belonging to its boundary dlZ 
are darker; 5 and S' are two superficial shells with thickness 
I outside 1Z; X — S U <S' U 1Z is the extended region; X is the 
exterior of X. 



Section VI contains a simpler proof for the area law (|3]) 
without assuming (P), but assuming ^ for all the range 
of e. In Section VII the above results for the ground state 
are generalized to other low-energy states (not necessarily 
eigenstates). Section VIII contains the conclusions. 



III. LOCALITY AND ENERGY 
FLUCTUATIONS 

A. Local interactions 

Consider a system with one particle at each site x £ 
£ of a finite s-dimensional cubic lattice L C Z s . The 
distance between two sites x, y £ C is defined by 



d(x, y) = max \xi 

l<i<s 



Vi 



(4) 



In the case of periodic boundary conditions or hybrids, 
this distance has to be modified with the appropriate 
identification of sites. Each particle x £ C has associated 
a Hilbert space with finite dimension q. 



The hamiltonian of the system can be written as 



(5) 



xEC 



where each term K x can have nontrivial support on first 
neighbors (y £ £ such that d(y,x) < 1). There is a 
constant J which bounds the operator norm of all terms 
\\K X \\ < J- (The operator norm of a matrix is equal to 
its largest singular value.) Translational symmetry is not 
assumed, so each term K x is arbitrary. The eigenstates 
and eigenvalues of H are denoted by 



H\* n ) = E n \* n ) 



The hamiltonian H satisfies the premises for the Lieb- 
Robinson Bound [U [22[ ■ Let X, Y be two opera- 
tors acting respectively on the regions X,y C £, with 
ll-^ll; ll^ll < 1- The distance between two regions is de- 
fined by 



d(X,y) = min d(x,y) . 



(7) 



The time-evolution of an operator in the Heisenberg pic- 



ture is X{t) 
states that 



Xe 



-iHt 



The Lieb-Robinson Bound 



\[X(t),Y}\\ <2\X\ 



(vt) 



[d(x,y)/2\ 



(8) 



where v = 2J5 S . When vt <§C d(X,y) the two operators 
almost commute. In other words, the dynamics gener- 
ated by H does not allow for the propagation of signals 
at speed much larger than v. A simple proof of the Lieb- 
Robinson bound ([5]) is provided in Appendix C. 



C. Average for the energy fluctuations 

For any region X C C and any integer I > 5 define the 
exterior, the boundary and the superficial shell as 

X = C\X = {x£C:x£X}, (9) 
dX = {x £ X : d(x, X) = 1} , (10) 
S = {x £ X : d(x,X) < 1} , (11) 

respectively (see FIG. 1). The hamiltonian Hx is defined 
as the sum of all terms K x whose support is fully con- 
tained in X. The eigenstates and eigenvalues of Hx are 
denoted by 



HxW„ 



(12) 



where the index n — 0, 1, . . . labels the eigenvalues in 
increasing order e„ < e n+ \. The sum of all terms K x 
which simultaneously act on both, X and X , is denoted 
Hi= H - H x - Hx, and has norm \\H X \\ < J3 S |<9A?|. 
(6) The expectation of any operator X with the ground state 
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is denoted by (X) = (* |-X"|*o)- Without loss of gener- 
ality it can be assumed that each K x is positive semi- 
definite, which implies 

(H x ) + (H x ) 

< {Hx+Hi + Hz) 

< tt[(H x +H X + Hx)WMq\ <8 tr*|tf X*o|)] 

< e + J3 s \dX\ + (H x ) , 



and 

eo < (H x ) < e + J T \dX\ . 
This can be sumarized as follows. 



(13) 



The energy frustration of the global ground state 
in a region X is, at most, proportional to the boundary 
dX. 



D. Observation of energy fluctuations 

For any value of e cu t define the operator 
dt 



2tt 



e i(E -u)t e iH x t e -iHt 



(14) 



where a = 10 4: v 2 /l. The action of Q onto the global 
ground state |^o) implements an approximate projection 
onto the subspace with energy lower than e cu t inside the 
region X, 



QI*o> = 



E 



\tpn)(tpn\ 



(15) 



This integral is the error function, which is a soft step 
function. In the limit where the softness parameter a 
tends to zero, the operator inside the square brackets be- 
comes a projector. The operator Q has non-trivial sup- 
port on the whole lattice C, but remarkably, it can be 
approximated by the operator 







e C ut 

dw 



dt 
2^ 



e i(E -ui) e iH s t e -iHj 



(16) 



which has non-trivial support only in the region X L) S. 
More quantitatively, the bound 



10-011 < l<*fe- 



(17) 



is proven in Lemma 1 (Appendix), using techniques sim- 
ilar to the ones in [l3|> [22|, [23[ . The fact that Q s» Q is 
solely a consequence of the locality of interactions and 
can be understood as follows. According to the Lieb- 
Robinson bound if i < l/v, any operator Y with 
support on X\S evolves to an operator Y{t) with ap- 
proximate support on X. Then e~ lHxt Y(t) e lHxt s» Y, 



or in other words, the unitary e lHxt e 



iHt 



in (0 ap- 



proximately acts like the identity inside X\S, or in other 



words e iH xt e -iHt py e tH s t e -iHx uS t^ ^yjjjgjj justifies the 
definition (fH 



The right-hand side of (TTBl is an average of unitaries, 
therefore ||0|| < 1. Then, the operators \Q\ and (I — \Q\) 
define a two-outcome generalized measurement on XUS, 
which tells whether the energy inside X is below or above 
e C utj approximately. 

Everything shown in this section for the ground state 
generalizes to all eigenstates. The action of Q onto |^ n ) 
is 



0|*n) = 



E 



e 2 " 

duj \lpn)\^n\ 



l*n) , (18) 



where e' cut = e cut +E n — Eo- Summarizing, for each eigen- 
state l^n) there is an operator Q which approximately 
projects onto the subspace with energy e n < e cu t inside 
the region X, by only acting on the exterior and the shell 
XUS. The degree of approximation increases with I, the 
width of S. The larger I is, the closer Q and Q are, and 
the smaller the softness parameter a is. 

The energy fluctuations of an eigenstate \*5f n ) inside a 
region X can be observed by measuring the exterior and 
a superficial shell inside the region, that is X U S (see 
FIG. 1). 



IV. SUPPORT OF THE GROUND STATE 
INSIDE A REGION 

A. Decay of correlations 

It is usually the case that, when the system is in the 
ground state, the correlation between two observables 
acting on different sites decrease with the distance be- 
tween the sites. Let T be a function which upper-bounds 
the connected correlation function of any pair of op- 
erators X, Y acting respectively on the disjoint regions 
X,yc£, with \X\ < \y\ and ||X||, ||F|| < 1, 



\(XY)-(X)(Y)\ <r(d(x,y),\x\) 



(19) 



For the argument of this paper, both, the decay with the 
distance d(X,y) and the scaling with size of the support 
of the operators \X\, are relevant. It is shown in [22| that 
any hamiltonian H with an energy gap above the ground 
state A = Ei — E > has 



r(Z,|*|)= Cl |*| 



-in 



with correlation length £ = lOv/A. To prove the area 
law for the entropy the following condition is needed. 

Assumption 1 The correlation functions for the ground 
state decay at least as 



r(i,|*|) = 



Cl 



(l-£ln\X\y 



(21) 
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where ci,£ and v > s are constants. 

Note that both, (j2"0)l and ([2~Tj) , have the same relative 
scaling of / and \X\, but assumption (|2ip is weaker than 
([20|) . Although the decay (l2~Tj) is polynomial in I, it is not 
the correlation function of a critical hamiltonian, where 
one expects T ~ Unfortunately, the argu- 

ment of this paper does not give an area law with such 
scaling in \X\. 

B. Energy fluctuations inside a region cannot be 
observed in its bulk 

For any region TZ C C and any integer I > 5 define the 
extended region as 

X = {x G C : d(x,K) < 21} , (22) 

which redefines ©, fllDJ and (JITJ (see FIG. 1). The 
region TZ can be considered the bulk of X. 

Suppose the existence of an operator Z with support 
in TZ such that 

Z|* ) « E l^nXV»n|*0> ■ 
n: e„<e cut 

This operator acts onto the ground state in a similar way 
as Q does, then the two operators are correlated 

(ZQ) a (Z) a (Q) , 

and their corresponding supports are separated by a dis- 
tance I. For the right choice of e cu t and large enough I 
the existence of Z is in contradiction with Assumption 1, 
therefore 

The energy fluctuations of the global ground state in- 
side a region X cannot be observed in the bulk of the 
region, that is TZ. 

In the following subsection, a quantitative example of 
this fact is given. 

C. Characterization of the support 

In what follows, the assignation 

e cut = 2J3 s |cU'| + e + 20v (23) 

is assumed in the definitions of Q and Q (|14I16[) . 

Definition of P For each eigenstate \tp n ) of Hx with 
e n < e cut + 2Qv consider the Schmidt decomposition [l[ 
IVVi) = YlilKiWn) ® \@n) with respect to the partition 
\a l n ) € Tin and |^) G Hx\n- Define P as the projector 
onto the subspace of Tin generated by all vectors \a l n ) 
defined above, symbolically 

P = supp TC {|?A„) : e„ < 2J3 s \dX\ + e + 40w} . (24) 



Let P = I — P be the projector onto the complementary 
subspace. Lemma 3 (Appendix) shows that the assigna- 
tion (|2U)) implies 

(Q) > \-m 3 e- 1 , (25) 
(P X Q) < 2\X\ 3 e~ l . (26) 

Recalling that the respective supports of P 1 - and Q are 
separated by a distance 2, one can invoke the decay of 
correlations (1191) without specifying the function T, 

(P X )(Q) - <P X Q) <r(Z,|ft|) . (27) 

The combination of p5]>. and (gT)) gives 

(P) > l-4T(l,\K\) , (28) 

for sufficiently large I, where 1/2 > r(/, |^|) > 6| A"| 3 e _z 
holds. Concluding, the support of the global ground state 
inside TZ is contained in the subspace characterized by P, 
up to some small weight (|28p . 

D. A renormalization group scheme 

The projector P defined above allows for certifiably- 
generating a low-energy effective theory for H : the hamil- 
tonian terms K x inside TZ can be renormalizcd as 

K x —— i PK X P . (29) 

The whole lattice can be divided in similar regions, and 
the transformation (|29l) performed in each of them. The 
fidelity between the effective and the original ground- 
states can be bounded with (|28l) . and increased by en- 
larging I. As explained in Section VI, one can also obtain 
arbitrarily-good fidelities for any low-energy state. 

V. ENTANGLEMENT IN THE GROUND STATE 

A. Energy spectrum 

In the previous section, a subspace which approxi- 
mately contains the support of the ground state inside 
a region has been characterized. In order to bound its 
dimension, an additional assumption is needed: if the 
boundary conditions of the hamiltonian are left open, 
the number of eigenstates with vanishing energy-density 
must not be exponential in the volume. 

Assumption 2 There are constants C2, t, 7, 77 such that, 
for any region X and energy 

e = 2J3 s \dX\+e + 40v , (30) 

the number of eigenvalues of Hx lower than e satisfies 

^x(e) = max{n : e n < e} 

< c 2 {T\X\)^ e - e ^+i\ dx \ . (31) 
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The area law is nontrivial when applied to regions TZ 
such that \dK\ < \TZ\, or equivalently \dX\ < \X\. 
In this case, the eigenstates with energy proportional 
to the boundary \dX\ (|30[) have vanishing energy den- 
sity e n /\X\. According to [23[, Assumption 2 holds for 
many systems that have an energy gap above the ground 
state. There are known hamiltonians which violate As- 
sumption 2 and have a gap, but when the boundary 
conditions are opened there appears a degeneracy for 
the ground state which is exponential in the volume 
23]. Massive free bosons and fermions satisfy Assump- 
tion 2. Contrary, massless free frermions violate it as 
Q ~ exp^e-eo)!*! 1 /*. 



The factor (t\X\) 



-y(e-eo) 



in (pJTj) can be understood 



with the following example. Consider the hamiltonian 



H x = 



1 




where the subindex x specifies in which site the matrix 
acts. The energy e£ {0,1,... \X\} counts the number of 
local excitations, hence the degeneracy is the binomial of 
\X\ over e, which can be upper-bounded by \X\ e . The 
constant factor (r|#|)''l a *l in (|3~Tj) is introduced because 
some hamiltonians with open boundary conditions have 
a degeneracy (or approximate degeneracy) which is ex- 
ponential in the size of the boundary. 

Consider again the Schmidt decomposition of each 
eigenstate \ip n ) with respect to the partition TZ and X\TZ 
(Definition of P). The dimension of the Hilbert space 
Hxvr. is q\ x \ n \ therefore the support of each \ip n ) on TZ 
has at most dimension q\ x \^. This and Assumption 2 
provide a bound for the rank of the projector P 

rankP < q \ x \^\ c 2 (t\X\)^ x ^ 2J3S+ ^+^ . (32) 



B. Entropy of an arbitrary region 

Consider a region 1Z C C being a completely arbi- 
trary subset of the lattice. It not need to be convex, 
full-dimensional nor connected. For any site x 

\{yeC:d(y,x)<2l}\ < (5l) s , 
\{ye£:d(y,x) = 2l}\ < 2s(5l) s - 1 , 



which imply 



1*1 < \K\(5iy , 

\dX\ < \dK\2si5iy- 1 
\X\TZ\ < \dlZ\(5l) s . 



(33) 



Let p-jz = tr c \7j|\I'o)(^'o| be the reduction of the ground 
state in 7Z, and Ai > A2 > • • • its eigenvalues in de- 
creasing order. Assumptions 1 and 2 imply (f2"Tj). (|2"8"|) 
and p2p. which impose the following constraints on the 



eigenvalues: for any integer I > 5, 

0(0 

fc=i 

4ci 



(34) 



9(1) = 



ln6(0 = |9^|2s(50 s_1 (72J3 ;5 + 7/)ln[T|^|(5/) s ] 
+ |9^|(50 s lng + O(ln|^|) . 

Now one can find the probability distribution Afe which 
maximizes the entropy (— ^AfclnAfc) given the above 
constraints. This is done in Appendix B with the follow- 
ing result. 

Result 1 The entropy of the reduction of the ground 
state inside an arbitrary region 1Z satisfies 



S(pn) < |0ft|(lO£ln|ft|) 

+ olidnimiiiy- 1 



ru + In q 



(35) 



C. Entropy of a cubic region 



Consider the case where the chosen region is a hyper- 
cube TZ — {x € C : < Xi < L}. One can proceed as 
before, but the bounds analogous to (l33|) are smaller, im- 
plying a smaller bound for the entropy. All this is worked 
out in Appendix B. 

Result 2 The entropy of the reduction of the ground 
state inside an cubic region TZ satisfies 

S(pn) < \dTZ\\xY\TZ\{ 1 2Ji s +7 1 +^\iYq)+0{\dTZ\) . (36) 



It is expected that the entropy of any full-dimensional 
convex region TZ obeys the same scaling (131)1) . 



VI. SIMPLER PROOF FOR THE AREA LAW 

An area law can be easily proven without Assump- 
tion 1, if Assumption 2 is extended to all values of the 
energy e, not only the ones satisfying (|3"0)) . Let TZ be the 
region where the entropy is estimated, and H-ji the sum of 
all terms of the hamiltonian ((5]) which are fully contained 
in TZ. Following the conventions of this paper, the eigen- 
states and eigenvalues are denoted by H-jz\ip n ) = e n \ip n ), 
where eo < e± < • • •. The strong version of Assumption 2 
tells that all the eigenvalues e„ satisfy 



n < c 2 {T\n\y {en - (io ^ an \ . 

The global ground state can be written as 



(37) 



(38) 



6 



where the coefficients fi n are non-negative and add up to 
one. It is shown in [l| that the entanglement entropy of 
| "Jo) is upper-bounded by the entropy of the /i-coefficients 

S(pn) < - ^/^nln/^n . (39) 

■n 

Locality implies equation (|13p. which can be written as 
J2^e n < e + J¥\dK\ . (40) 

n 

Maximizing the right-hand side of (|39D over the proba- 
bility distribution jj, n and the numbers e„ subjected to 
the constrains (I3T1) and (|4"0")) gives 

S(pn) <2J3 s 7 |9^|ln|^| , (41) 

the area law. This calculation is made in Appendix D. 

VII. ENTANGLEMENT IN EXCITED STATES 

The entanglement properties of excited states have also 
been studied. In references |24T - |26| the motivation was to 
study the robustness of the area law for the entropy of 
black holes. They show that in systems of free bosons, 
the entropy of some low-energy excited states scales at 
most like the area. In [27| the entropy scaling in one- 
dimensional spin systems is analyzed. They show that 
some excited states have entropy proportional to the vol- 
ume, but low-energy states obey an area law. All this 
work is for integrable systems. In what follows, we ad- 
dress the general case. 

Sometimes, low-lying excited states \^ n } have corre- 
lation functions similar to the ones of the ground state. 
The single-mode ansatz for excitations with momentum 
k is 

|#n cx 5> fa - fc £ x |# ) > ( 42 ) 

X 

where Z x is an operator acting on site x such that (Z x ) = 
0. If X, Y have support on finite regions and the volume 
of the system \C\ tends to infinite, then the correlation 
function ([T9"f for the state (l4"2"j) is the same as for \^o). 
The same happens to excited states containing a small 
number of single-mode excitations. Examples of single- 
mode excitations are: spin waves, free bosons and free 
fermions. In this section it is shown that such excited 
states obey an area law similar to the one for the ground 
state. Actually, this is done with a bit more generality. 

Consider an arbitrary superposition of eigenstates with 
bounded energy 

n:E n <E Ia 



In this case, the correct assignation for e cu t in the defi- 
nitions of Q, Q and P (H [H EH) is 

e cut = 2 J3 S \dX\ + e + 20v + E m - E . (44) 

Applying Assumption 1 to the state (j4"3"]l . the arguments 
follow exactly as for the ground state. Repeating the cal- 
culation of the entropy for a cubic region 1Z, and keeping 
track of the term proportional to (E m — Eq) one obtains 

S(tr ATC |$X$|) < |dft|ln|ft|( 7 2J3 s + 77 + 4£mg) 
+ (E m ^E )(\nT\U\) 1 - ^ Cl ^ +3 
+ <D{\dK\) . (45) 

VIII. CONCLUSIONS 

It is shown that ground states and low-energy states 
obey an area law for the entropy, provided two conditions 
hold: (i) the state has a sufficient decay of correlations, 
and (ii) the number of eigenstates with vanishing energy- 
density is not exponential in the volume of the system. 

A universal property for local hamiltonians is also here 
established. The energy fluctuations of eigenstates inside 
an arbitrary region can be observed by measuring the 
exterior and a superficial shell of the region. This extends 
to any pure state that can be written as a superposition 
of eigenstates with similar energy. 

Some thermodynamic quantities at finite temperature 
only depend on the density of states. Examples are: free 
energy, (global) entropy, heat capacity, etc. This pa- 
per establishes a relation between these thermodynamic 
quantities and ground-state entanglement. 
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Appendix A: proofs 

Lemma 1. Let Q,Q be the operators defined in (fT4"]) . 
(dU), then 



\Q-Q\\ < x\ e 



3 „-l 



(Al) 



Proof First, express Q and Q with a single integral, by 
using the identity 



Z7T 



dt e 2 



2?r 0+ - it 

Second, define the operators 

Ho = Hgugr — Hs — H$> 
H\ = H — Hx — H x , 



which respectively act on the regions Ho, Hi C C. Note 
that dCHo,-Hi) =1-4, \U \ < \X\ , \Hi\ < 5 S \X\, and 

= II e i{H-H-,)t e -iHt _ gi(H—Hx—Ho)tg—i(H—Ho)t || 

This, the triangular inequality, Lemma 2, and the Lieb- 
Robinson bound ©, give 



\Q-Q\ 



< 2 dt 



1 
2~7rf 



\ e iH x t e -iHt _ e iH s t e -iH xus t\ 



2 / dt — 2 

to 27rf 

3 t2cs pto 



2\X\ A J Z 5 S f l ° 1 / , / - yvv! 
/ dt - I dt 2 I dh 



7T 

4 



: e 2 







t 2 ( wtl U'/2-2j 



L//2-2J! 



< 



toy/a 
\X\ 3 (^ ) Li/2j 



7r5*(Z-l) [l/2\\ toVo" 



e 2 



Putting to — [I /2J I \e 3 v) and using Stirling's approxima- 
tion 



(vt ) 



L//2J 



< e 



[1/2] In 



Tm] < e 



1-/ 



LZ/2J! 

Putting a = 10 4 v 2 /l > 2lt^ 2 one obtains ([AT]) . □ 

Lemma 2. Let H,X,Y be hermitian matrices and 
t > 0. then 



(A2) 



\\ e i{B-X)t e -iHt _ e i(H-X-Y)t e -i(H-Y)t\ 
< [dt 2 f dt! || [Xfa), Y] || , 

Jo Jo 
where X(t) = e iHt X e - lHt . 

Proof If f(t) is a diffcrcntiable function with /(0) = 
then f(t) = f Q dtif'(ti). This implies the following two 
equalities. The following two inequalities are a conse- 
quence of the triangular inequality for the operator norm. 

i(H-X-Y)t e -i(H-Y)t e Wt e -i(H-X)t _ jj|| 
H-X-Y)l 2 r _ • v c -i(H-Y)t 2 JHt 2 



- e -iHt2 e i(H-Y)t 2 \ 













If 










H 


- e~ l 








fdt 2 


" J 











fdt 2 


= J 











fdt 2 


~ J 






dt x e'< H ' Y ^ ^,X(t x )\e< H - Y ^ 
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Lemma 3. The operator Q defined in ([T5| with e cu t = 
2J3 5 |<9A'| + eo + 20v, and the projector P 1 - defined by 
(|2l]). satisfy 



(Q) > --2|Af e -* , 



(P X Q) < 2|#| d e 



(A3) 
(A4) 



Proof . The positive operator 



M = I duj (2™y 1/2 \i> n ){ipn\ (A5) 

(A6) 
(A7) 

(A8) 



allows for writing equality (|15[) as 

Q|*o) =M|* ) • 

The two projectors 

n: e„ <e cu t ±5 

with (5 = 20i>, satisfy 

M_ - e~~*I < M < M+ + e _, I , 

where we have used that < e~ l . The positivity of 

M and the second inequality in (|A8[) imply 

M 2 < (1 + 2e~ l )M + + e~ 21 . (A9) 

A worst-case estimation gives 

(H x ) > (M_)e + (I - M_)(e cut - «J) . (A10) 

Performing the assignation e cu t = 2J3 S |9A'| + cq + 5 in 
(IAT0|) and using bound (US]) one obtains (M_) > 1/2. 
The combinations of (fI7|). (|A"6l) and (|A"8]l gives (1A"3)) . 

Using Lemma 1 and (|A6I) . the Cauchy-Schwarz in- 
equality, bound (IA9|) . and the definition of M + and P- 1 , 
one obtains respectively the following chain of inequali- 
ties: 



3 — £ 



(P X Q) < (P x M) + |#| d e 
: (P ± ) 1 / 2 (P ± M 2 P ± ) 1 / 2 + \X 



-lyplir d_L\ i „-2il 1 / 2 , I v \3„-l 



< [{l + 2e- l )(P ± M+P 



r'+\x\*e- 



< 2\X\ 3 e- 1 



which is [KM. 



(All) 

□ 



Pk 
Pk 



Appendix B: calculation of the entropy 

1. Entropy of an arbitrary region 

Consider the probability distribution defined by 
1 - 0(h) 



e(z ) 

6(1 - 1) - 6(1) 



0(0-9(1-1) 



for 1 < k < Q(l ) , (Bl) 
for Q(l - 1) + 1 < k < 6(0 , 



for every integer I > Iq = 2£ In \1Z\, and 

«" " g-W ■ <B2) 

ln©(0 = |a^|2s(50 s_1 (72J3 s + r;)ln[T|^|(50 ;i ] 
+ |a^|(5/) s lng + C(ln|^|) . 

This distribution is uniform in blocks of the maximum 
size that constraints ([M]) allow. Then, it is the distribu- 
tion satisfying (|34p with maximum entropy. The upper- 
bound on the entropy of pk gets simplified by using the 
substitutions 6(0 - 9(Z - 1) < 6(0 and 



6(1-1) -6(1) < 



2^+3 



(i - i\n\n\y+ 1 

Using this, one obtains 



(B3) 



+ olidnKHniy- 1 ] 



2. Entropy of a cubic region 



(jJ3 s + rf) +lnq 



(B4) 



Consider the case where the chosen region is an hyper- 
cube 1Z = {x £ C : < Xi < L}. It is easy to calculate 

\n\ = l s , 

|072.| = 2sL s - 1 . 

Following definitions ([2"2l ITU1) one obtains 

\X\ = \H\(l + Al/L) s , 
\dX\ = \dK\(l + 'll/Ly~ 1 , 
\X\1Z\ < \dX\2l . 



Consider the probability distribution (|B1[) with 6(1) given 
in (JB2J but 6(0 defined as 

ln6(0 = |c>^|(l + 4Vi) s_1 (72J3 ;s + 77)ln[r|^|(l+4Z/L) i 
+ |^|(l + 40^) s_1 2nn g + O(ln|^|) . 

Using the same tricks as above one obtains the following 
upper-bound for the entropy of pk, 

-^Pklnpk < \dK\ln\K\(j2J¥+? 1 +^lnq)+0(\dll\) . 



Appendix C: The Lieb- Robinson bound 

Let X, Y be two operators with support on the re- 
gions X, y respectively, and L = d(X,y). Let Z be an 
arbitrary operator and F(t) = [X(t),Z], where X(t) = 
e lHt Xe~ lHt and H is the hamiltonian ((SJ). Using the 
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Jacobi identity [[X,Y],Z] + [[Y,Z),X] + [[Z,X],Y] = 
twice one obtains 

d t F(t) = i[[H,X(t)],Z] (CI) 
= -i[F(t),H]-i[[Z,j: xeZ K x ],X(t)} 
= i[A,F(t)]+tZ aez [Z,[X{t),K x ]] , 

where Z = {x : [K x , Z] ^ 0} and A = J2 xe c\z K *- The 
above is equivalent to 



St (e^ At F(t) e lAt ) = i J2 e~ lAt [Z, [X(t),K x 



J. At 



x£Z 



which can be integrated 

e - tAt F(t) e lAt 



(C2) 



= F(0) + z^ [dhe-^iZ^XihlK^e^ . 
xez Jo 

The triangular inequality for the operator norm gives 

\\[X(t),Z}\\ (C3) 
< \\[X(0),Z}\\+2\\Z\\J2 [dh\\[X( tl ),K x }]\\ . 

Define g x (t) = \\[X(t), K x }\\ and use |C3]) with Z = K x 
to obtain 

g x (t) < g x (0) + 2J f dti g x >{t x ) . 

X ':d(x,x')<2 

If r = d(x, X) > 2 then 5^(0) = 0. The above can be 
iterated n = [(r — 1)/2J times 

pt pt n pt 2 
g x (t) < v n max / dt n / dt n ^ 1 ■■■ / dt x g x >(h) , 

x'\ d( X , X ')<2n Jq J q J q 

where v = 2J5 S and \{x' : d(x,x') < 2}| = 5 s . For any 
x' the bound g x >(0) < 2J\\X\\ holds, then 



g x (t) < 2J\\x\ 



( ^L(r-l)/2J 

L(r-1)/2J! 



(C4) 



Note that |{a; : [K X ,X] ^ 0}| < 5 S \X\. This and the 
bound (|C4| can be substituted in (|C3j) with 2 = Y, 
giving 



[X(t),y]||<2|A-|||X||||y| 



This is the Lieb-Robinson bound. 



[L/2\\ 



Appendix D: Simpler proof for the area law 

Let us obtain an upper-bound for the entropy 



S = - /i w In /i r , 



(Dl) 



of any probability distribution /i„ subjected to the con- 
straints 

n < c 2 (r|^|) 7(e "- eo)+ " |aK| , (D2) 
XXe„ < e + jy\dTZ\ . (D3) 



The original problem has a finite range for n, but relax- 
ing this fact cannot diminish the largest entropy (|Dlj) . 
Isolating e n from (|D2|) gives 



ln(n/c 2 ) ??. . 
e n > , , „ -~\n\+e , 
7ln(r|^|) 7 

which substituted in (|D3|) gives 
^/x n lnn < |dft|ln(r|ft|) (J3 s 7 +77)+lnc 2 = c 3 . (D4) 



The new constant C3 is define through this equality. For 
any a > 0, the entropy (|Dlj) can be upper bounded as 



S < max -) jx n lnpi n + a(C3- }^ fi n lni 



(D5) 



This upper-bound holds because the second term be- 
tween the square brackets is always positive. In what 
follows, we perform the maximization (|D5|) with no con- 
strains over {/i n }: normalization, constrains (|D2[) and 
(|D3I) are ignored. This relaxation cannot decrease the 
right-hand side of (|D5[) . Deriving the right-hand side of 
(ID5[) with respect to /i„ and equating to zero gives 

— In /i n — 1 — a In n = , 

which implies /i„ = rT a je. Substituting this into (|D5[) 
gives 



5 < 5>-/ 



e + ac 3 



Choosing a = 2 and replacing C3 gives 
5 < 2J3 s 7|97e|ln(r|^|) +2r]\&R\ +21nc 2 + 7r 2 /(6e) 



